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Abstract 

A spin angular gradient approximation for the exchange correlation magnetic field in the density functional formalism is 
proposed. The usage of such corrections leads to a consistent spin dynamical approach beyond the local approximation. The 
proposed technique does not contain any approximations for the form of potential and can be used in modern full potential band 
structure methods. The obtained results indicate that the direct 'potential' exchange in 3d magnets is rather small compared 
to the indirect 'kinetic' exchange, thus justifies the dynamical aspect of the local density approximation in 3d metals. 



In density functional theory (DFT) two magnetic terms, responsible for any spin ordering, coexist on both in- 
traatomic and interatomic length scales and produce a very successful description of static properties in 3d magnets. 
One of these terms is the static non-local spin current, which is associated with the kinetic energy ('kinetic' field 
(KF)), whereas another one is the large exchange-correlation magnetic field (EXF). While enjoying great success in 
weakly correlated materials, the standard local spin density approach (LSDA) for the exchange correlation energy 
is believed not appropriate for strongly correlated materials. However, the standard assumption of the LSDA - a 
proportionality of EXF and local magnetization m(r) is, in turn appropriate only for strongly correlated systems 
(narrow bands), when the magnetic effects related to the 'kinetic' exchange can be omitted. The consequence of 
this local approximation is that the large EXF is spin dynamically inactive in the LSDAEl. This feature of LSDA 
restrains the possible range of directions for EXF, which immediately affects the spin-spin correlation function and, 
in turn a whole range of dynamic and thermal properties. Hence, it is highly desirable that both sources of internal 
magnetic field are treated with no constraints on their directions. A general non-local approach should be used to 
treat adequately the dynamic aspect of the interrelation between 'kinetic' (indirect) and 'potential' (direct) magnetic 
sources. The influence of non locality of EXF on spin dynamics (SD) of magnets is essentially unknown and was not 
studied in real materials. In this paper we undertake an attempt to overcame this problem and propose a spin angular 
gradient approximation (SAGA) for EXF and apply it to the calculations of spin wave (SW) stiffness in 3d magnetic 
metals. 

The evolution of the many-electron inhomogeneous system in the presence of the external field is defined in a 
unique way by the time-dependent one-particle density matrix p a[3 (r, r,i) = ^-0/3 ( r ji)V'a ( r ^)^) ; where %p a (r, t) is 
the annihilation operator for the electron at the point r with spin projection a at the instant time t. Equivalently, 
one can introduce the charge n (r,t) — Trp (r, r, t) and magnetization m (r,t) = Trp (r, r, t) cr densities, where er are 
Pauli matrices. Starting from the Schroedinger equation for the many-electron system one can formally obtain the 
exact set of equations for these quantities: 

m(r,t) = 7m(r,t) x B ext (r,t) + -Vl T ,{p a p (r, r' ,t) a 0a - c.c.) r < =r , (1) 



n(r,t) = l -V 2 rr ,( Paa (r, r',t) - c.c.) r , =r (2) 

where B e2:t (r,i) is the external magnetic field. These equations are not closed since, generally speaking, the quantity 
V r 'yO Q/3 lr, r',t) r ,_ r cannot be described directly in terms of p af3 (r, r,t) (see the discussion of the kinetic energy term 
in DFTtl). Introducing the wave functions of Kohn-Sham quasiparticles (p l/a and making an adiabatic approximation 
which assumes they evolve rapidly in comparison to the spin degrees of freedom, we can obtain a closed set of equations 
describing charge and spin dynamics: 

m(r,t) = 7m(r,t) x B tot (r,i) + (3) 

OCC 

occ 

= ^ V r(H <P*ua M V r^ Q (r,t) - c.c), (4) 



-V 2 + V ex t + Vh) 5 a p — BextCaP + f 
2 / °P/3a 



Vvf) (r,i) (5) 



where the sum is over the occupied Kohn-Sham states; V ext and Vh are the external and Hartree potentials, re- 
spectively, E xc is the exchange-correlation energy, B toi = B e:Et + TS XC is the total magnetic field including EXF 
B xc = — 5E XC /Sm. The Kohn-Sham wave functions and the corresponding energies e v (t) depend on time due to the 
time-dependence of the densities and external field (the latter is supposed to be weak relative to the characteristic 
electron energies). 

The second term in Eq.(||) can be conveniently cast in the following torque form 

- = 7m(r,t) x [B fem (r) + B xc (r) + B ext (r)] (6) 
= 7m(r,£) x B tot (r) 
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where the total field acting on the electronic spin at point r is 
r> < v SE ST 6E XC 

B tat r = - — = - - - +Bext{r), 7 

om(r) om(r) om(r) 

with Bfcj ra (r) = VpVcr/mp, p = J] (fttfiui d/cft = <9/9i + vV and v is a velocity. From Eq.(^J) directly follows that 
mdm/dt = 0. Below we use B ext = 0. 

The basic assumption of LSDA is that E xc is obtained for the homogeneous electron gas (HEG) model for a given 
magnetization density m, and assumed to apply to real external potentials. In this local approximation H xc (r) ~ 
m(r), hence assuming very strong EXF compared to KF and smallness of dynamics associated with EXF. Let us now 
eliminate this requirement and generalize a shape of EXF for the spiral magnetic configuration (spin-density wave 
(SDW)) with the amplitude m and the wave vector Q. Hence, the term m(r,i) x B xc (r) will contribute to Eq.(||). 

The SDW state is characterized by anomalous averages s p — ^ c p + Q/2t c p-Q/ 2 i) wnere c po-> c p<x are the oper- 
ators of creation and annihilation of electrons with momentum p and spin rarojection a. Similar to the Gorkov- 
Nambu formalism in the theory of superconductivity (see, for instance, Ref.u) we introduce the spinor operator 
tp p = ( c p + Q/2f c p-Q/2l)- Then the Hamiltonian can be expressed as 

H = £ tf p Mp ^EE^W (</Wp) (V£- q V) (8) 
p q#o pp' 

where v c (q) = 47re 2 /q 2 V, V is a volume and h p — 9 p + r p a z — A p <r x . Here 

9 p = \ ( £ p+Q/2 + £ P -Q/ 2 ) = P 2 /2 + Q 2 /8 ~ M> (9) 
t p = \ ( e P+Q/2 - ^p-Q/2) = PQ/ 2 

where e p = p 2 /2 — [i is the energy of the free electron and 2A P is the antiferromagnetic gap, related to the formation 
of SDW. In the Hartree-Fock approximation (HFA) the latter is determined through the relation 

a p = I> c (p-p>p'- (10) 

p' 

Now we will replace v c by the effective Stoner exchange splitting I — {V/) xc — V} xc /) / (nj — njj , where V° xc — 
d irie exc ) /drier- Then, Eq.[10] can be replaced by A = / (n-\ — njj /2, where A does not depend on p. 
The 'bare' Green function in the Matsubara representation has the form 

G{tw m ,p) = - — = 7 \ / !i -^r (11) 



ioj m - h p 



where £ p j ^ = P =F E p is a quasiparticle spectrum for SDW with E p = yT p + A 2 . From Eq.|Tl| one can find the 
occupation number matrix 

( T-nUz — Aa T \ ( tv-ov — Aav \ 
2N P = {l + -B-^ * j /pT + ^ _ _J^_ « j /pi (12) 

where / p(T = / (£ pcr ) is a Fermi function. Then for the Fock contribution to the exchange-correlation energy we will 
have 



Epock 
^Fock 



p (2) 
^Fock 



\ £ vc (P - P') Tr [JV(p)JV(p')] = 4i & + 41*, (13) 
pp' 



pp'er 



^p^p' 



^Mp-p')/«V,(.-^± 

PP' V lr 
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Expression ( |13[ ) is suitable for the calculation of the exchange energy in SDW state for any Q with the corresponding 
equation for the chemical potential /i 



n = TtJ2n p = J2U«- 



(14) 



Let us consider an important case of small Q. The expansion of Eq.[l4| gives a correction of order Q 2 for /i = [i — Q 2 
One has 



<*M = Mq - Mq=o 



Q 2 



8F(n T ,nJ 



(15) 



where F = (pf\ + Pf\ )I( n i, n-i )/27r 2 is a dimensionless Stoner enhancement factor, pp a = (6ir 2 n a ) 



1/3 



Expanding EqJl3 up to Q we have the following expression 



ock 



e 



(p F1 +P Fl )-Q 2 



12F 2 



(16) 



Now let us consider the correlation effects. In the Gell-Mann-Brueckner theoryi the correlations can be taken into 
account by using the polarization operator 



II(iw,q) = -Tr^T^G(p + q,i£ n + ^ n )G(p,i£„). 



(17) 



With this function the correlation part of Sl-potential can be expressed as: 

oo 

1 + v c (q) II (iu, q) 



n,. 



1 + v c (q) n Q=0 (iu, q) 



v c (q) [n (iu, q) - n Q=0 (iu, q)] 



where we wrote only Q— dependent part of the correlation energy. Substituting Eq^ into Eq.^we have 



n ( ^ q) = iWi + ^ 



A 2 



^""^ /per /p+qo- 



2 1 



p ' p+q 



iuJ + Cp+qcr - £p<x 
A 2 \ / pT - /p+qj. 



EpE p+ q 



iuj + £ 



p+ql 



with 



n Q=o (^, q ) = ^0 i + - 



2 / «*> fl 



In 



+ q 2 /2 + gp Fo 
+ <7 2 /2 - qpFo 



1 






2^ 


Pf<j ~ 


U v 



(18) 



(19) 



In 



- q 2 /2 + qp Fo 



iu - q 2 /2 - qpFo 
(20) 



where In means the main branch of the logarithm. 

Now the corresponding exchange-correlation addition to the spin wave spectrum at finite Q can be written as 



6u Q = — [E SD w(Q) - E SDW (0)] 
M 



(21) 



where Esdw(Q) is the total energy of the spin spiral and M is a magnetic moment of the cell. 

For the electronic gas of high density the main contribution into integral ( |l8| ) comes from the region of small q. 
At the same time from Eq.(|l^) one can see that the interband contribution contains additional factor q 2 compared 
to the intraband one. Hence, it seems reasonable to take into account terms of q 2 order for the correlation effects in 
intraband transitions. 

Below we consider non-relativistic case only, when exchange-correlation energy is invariant with respect to the 
rotation in spin space. Then, for the weakly nonuniform distribution of spin density, one can write 



dr{ne exc (n T ,n x ) + X(n h n l ) D} . 



(22) 
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where D — (V Q e^) (V a ep) — (V9) 2 + sin 2 9 (V</5) 2 is the rotational invariant of lowest order. Here e = m/ |m| = 
(sin cos ip, sin 9 sin ip, cos 9) . 

For the variation of E exc in the spiral structure with e(r) = (cos Qr, sin Qr, 0) one can write E® xc — E^. — VXQ 2 . 
Now, one needs to estimate the energy of SDW with an amplitude |m| = n-f —n± and Q — > 0. This problem has been 
considered by Herringtl in HFA. The last one, however, seems not to be useful in the calculation of the total energy of 
the SDW due to an instability at Q. = 2/cpH At the same time, when screening is taken into account this instability 
is destroyed (see, for instance, Ref.cl). Hence, it is essential to add screening effects to the Fock contribution to E exc . 
So, the following expression for the function of A in Eqj2^ can be obtained 

e 2 / 1 4\ e 2 
A (n T> nJ = [f~3J ( V ^ Pf ^ ~ V ^cPfi) + 7^5^ ■ (23) 

This expression properly takes into account the main effect - shift of chemical potential ( |l5| ) and exactly corresponds 
to Eq.|l^ in the Fock approximation. 



Let us discuss the consequences of the addition (23) to exchange-correlation energy. First of all, EXF B e3:c (r) = 
— 5E exc /5m(r) is now noncollinear with respect to the local magnetization, so, the variation of Eq.|22| leads to appear- 
ance of a new term 

B2 xc (r) = - (6 01 - e/3 e 7 ) V Q (AV a e0(r)) . (24) 
m 

The noncollinearity of B exc with e corresponds to the 'direct' exchange processes, which are absent in LSDA. This 
field can be directly included in any traditional full-potential LSDA technique. 

Another important consequence of this non-locality is the appearance of an additional contribution for SWS of 
fcrromagnons 



D=± 

M 



E SDW (Q) - E SDW (0) 



For this contribution we have 



(25) 



5D = ^ [ dr\(n hni ), (26) 

with integration over the whole elementary ccJl. This correction represents an addition to the value of D which was 
obtained in the random phase approximationp. 

The numerical calculations were first performed for HEG model. The calculation of the Q dependence of the 
correlation energy according to Eq.(16) revealed that this contribution is relatively small and does not have any 
peculiarities as a function of Q and has monotonic behavior. In our approximation it could indicate that this term 
contributes to both intraatomic and interatomic range scales in similar fashion. 

Further, we applied the technique described above for the studies of SD in ferromagnetic (FM) Fe and Ni. These 
systems are well studied and might be used as prototype systems for 3d magnetism research. Below we employ a 
non-collinear version of the full-potential linear augmented plane wave method (WIEN code) using both LSDA and 
generalized gradient correction approaches. In this approach the magnetization density is treated as a continuous 
vector quantity. The parameters of calculations (number of k-points, plane waves and etc) were chosen to provide 
a convergence of SWS about 3-5%. Unfortunately, we could not compute the energy from Eq.(16) for arbitrary Q 
with real DFT wave functions and our study of spin waves in Fe and Ni was performed only for small Q, i.e. for 
SWS from Eq.|5[ Our standard LSDA calculations (without SAGA) produced 239 meV/A 2 and 692 meV/A 2 for 
D in Fe and Ni respectively (see, also RefJl3). These results are in agreement with previous theoretical calculations 
and are close to the experimental numbers. The non-self-consistent addition of non-local correction from Eq.^6] does 
not change significantly these results, producing a positive correction of 13 meV/A 2 and 45 meV/A 2 for Fe and Ni 
respectively (so, the total D become 253 meV/A 2 and 735 meV/A 2 ). The addition of Eq.^2] for the self-consistency 
procedure and the total energy calculations modifies these results. The self-consistent account of SAGA increases the 
total D up to 271 meV/A 2 , with the increasing EXF contribution of 21 meV/A 2 . The effect of the increase of the 
total D was also found in Ni where EXF contributes 94 meV/A 2 to the total self-consistent value of 782 meV/A 2 . 
This result already indicates that SD in 3d metals exists mostly due to 'kinetic' (indirect) exchange with relatively 
small contribution from 'potential' (direct) exchange. Whereas the absolute value of the latter term is rather large, its 
non-local (dynamic) part is much weaker compared to 'kinetic' spin current dynamics. Overall the addition of SAGA 
leads to an increase of the total D from both the direct SAGA addition and the indirect modification of 'kinetic' 
contribution (effects of self-consistency). Whereas this result is obtained only for small Q, we expect that at larger Q 
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this effect will be more pronounced. Also one can argue that the addition of diagrams beyond the GW approximation 
can be important due to their influence on Stoner exchange splitting. Both LSDA and GW calculations produce much 
larger splitting than the one observed on experiment, hence one can expect that the other diagrams will correct it 
and as a result will generate a smaller SWS. 

In summary, we estimated the non local part of the exchange correlation magnetic field in real magnets. Our 
results for Fe and Ni indicated that in spite of the large absolute value of the exchange correlation magnetic field, the 
non-local part of this field, which contributes to microscopical spin dynamics is rather small compared to the 'kinetic' 
spin current term in 3d ferromagnetic metals. We expect that this correction can be important in the interstitials 
of weak magnets, systems with strong interaction between the local moment and conduction electrons, and magnets 
with non-collinear ordering. 
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